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ABSTRACT 

We produce mock angular catalogues from simulations with different initial power spec- 
tra to test methods that recover measures of clustering in three dimensions, such as the 
power spectrum, variance and higher order cumulants. We find that the statistical 
properties derived from the angular mock catalogues are in good agreement with the 
intrinsic clustering in the simulations. In particular, we concentrate on the detailed 
predictions for the shape of the power spectrum, P(k). We find that there is good 
evidence for a break in the galaxy P{k) at scales between 0.02 < k < 0.06 ftMpc -1 
using an inversion technique applied to the angular correlation function measured from 
the APM Galaxy Survey. For variants on the standard Cold Dark Matter model, a 
fit at the location of the break implies Slh = 0.45 ± 0.10, where f2 is the ratio of the 
total matter density to the critical density and Hubble's constant is parameterised as 
Hq = 100 hkm s _1 Mpc _1 . On slightly smaller, though still quasi-linear scales, there is a 
feature in the APM power spectrum where the local slope changes appreciably, with the 
best match to CDM models obtained for VLh ~ 0.2. Hence the location and narrowness 
of the break in the APM power spectrum combined with the rapid change in its slope 
on quasi-linear scales cannot be matched by any variant of CDM, including models 
that have a non-zero cosmological constant or a tilt to the slope of the primordial P(k). 
These results are independent of the overall normalization of the CDM models or any 
simple bias that exists between the galaxy and mass distributions. 
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1 INTRODUCTION 

Angular catalogues of galaxy positions provide us with pow- 
erful constraints on theories of structure formation in the 
universe. The APM Galaxy Survey covers 4300 square de- 
grees on the sky and contains over 2 million galaxies to a 
limiting apparent magnitude of b.j < 20.5 (Maddox et al 
1990a,b,c; 1996). The shape of the angular correlation func- 
tion measured from the survey at scales of 6 > 1° indicates 
that the universe contains more structure on large scales 
than is predicted by the standard Cold Dark Matter sce- 
nario (Maddox et al 1990c). 

Whilst this result is confirmed by the largest redshift 
surveys currently available {e.g. Efstathiou et al. 1990a, 
Saunders et al 1991, Vogeley et al 1992, Fisher et al 1993, 
Tadros & Efstathiou 1996), measurements of correlations in 
3D catalogues are still noisy on scales r > 10/i _1 Mpc. Only 
after the completion of the Sloan Digital Sky Survey (Gunn 



& Weinberg 1995) will a 3D catalogue contain the same or- 
der of magnitude of objects as the APM Galaxy Survey. 

An additional complication in redshift catalogues is that 
the pattern of galaxy clustering is distorted by the peculiar 
motions of galaxies (Kaiser 1987). This effect can boost the 
amplitude of the measured two-point correlations by any- 
thing between a factor of 1 — 2 on large scales depending 
upon the survey and the method of analysis (see Table 1 in 
Cole, Fisher and Weinberg 1995). 

Whilst the next generation of redshift surveys will un- 
doubtedly provide a wealth of new information that is not 
available in angular catalogues, it is important to take full 
advantage of the large number of galaxies and volume sur- 
veyed in the angular catalogues (such as the APM and the 
parent catalogue for the Sloan Survey when it is complete) to 
extract information about the correlations on large scales. 
Under certain assumptions, deprojection algorithms to re- 
cover the 3D correlations in real space have been developed 
for multi-point correlation functions (e.g. Groth & Peebles 
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Figure 1. The power spectrum measured for APM Survey galax- 
ies is shown by the open circles with the lcr scatter on the mean 
averaged over the survey split up into four zones. The solid line 
shows the linear power spectrum estimated from this as described 
in the text. The dashed line shows the power spectrum measured 
from the evolved simulation, APM2(a). 



1977, Fry & Peebles 1978, Peebles 1980), for J-order cumu- 
lants of counts in cells (Gaztanaga 1994, 1995; hereafter G94 
and G95) and for the power spectrum (Baugh & Efstathiou 
1993, BE93; 1994, BE94). In this paper we present tests of 
these algorithms by constructing angular catalogues, with 
the same selection function and angular mask as the APM 
catalogue, from large numerical simulations. We use sets of 
simulations that have been evolved to have a 3D power spec- 
trum that matches closely the APM form recovered by BE93 
and also simulations of CDM models. 

For our present purposes, we are concerned with testing 
for the presence of any systematic biases that arise from the 
projection process itself rather than from the actual con- 
struction of the APM Survey or corresponding angular cat- 
alogue (some of these problems are addressed in detail by 
Maddox et al 1996) 

The outline of the paper is as follows. In Section 2 we 
describe the N-body simulations used to make mock cata- 
logues in Section 3. We present and test the recovery meth- 
ods in sections 4 and 5. In Section 6 we discuss our results 
and present the conclusions. 



2 N-BODY REALISATIONS OF APM GALAXY 
CLUSTERING 

A discussion of the production of evolved N-body simula- 
tions that have the same power spectrum as that measured 
for APM Survey galaxies (BE93, BE94) is given in Baugh 
& Gaztanaga (1996; BG96). In this Section, we briefly sum- 
marize the approach taken and list the parameters of the 



Table 1. Simulation parameters. The third column gives the di- 
mension of the FFT mesh used in the long range force calculation. 
The last column gives the softening length of the gravitational 
force, e in h~ 1 kiloparsec units. 



run 


number 
of particles 


mesh 


(ft^Mpc) 


e 

(h^kpc) 


APM1 


160 3 


256 3 


440 


115 


APM2(a) 


200 3 


256 3 


600 


156 


APM3(a)-(e) 


126 3 


128 3 


400 


104 


SCDM2(a) 


200 3 


256 3 


600 


156 


SCDM3(a) 


126 3 


128 3 


400 


104 


LCDM3(a) 


126 3 


128 3 


400 


104 



N-body simulations that are used to make mock APM cat- 
alogues in Section |^. 

The first step is to estimate the linear power spectrum 
from the measured power spectrum of APM Survey galax- 
ies. This requires assumptions to be made about the cos- 
mological model and the form of the bias, if any, between 
fluctuations in the light and the mass distributions (Kaiser 
1984). In this paper we consider a spatially flat universe with 
the critical density fl — 1 and zero cosmological constant. 
We assume that there is no bias between light and mass, 
i.e. that light traces mass, for simplicity. The validity of this 
assumption is not important for the purposes of this paper, 
which are to generate a particular distribution of points in 
three dimensions and to determine how well the N-point cor- 
relations in 3D can be recovered from a projected catalogue. 
There is evidence that the relative bias between mass and 
light is small on large scales from the hierarchical scaling of 
higher order moments of galaxy counts in the APM Survey 
(G94), though this does not appear to be the case on smaller 
scales (BG96). 

The linear power spectrum is obtained from the evolved 
power spectrum using the transformation of Jain, Mo & 
White (1995). This transformation is based upon a sugges- 
tion by Hamilton etal. (1991) that a universal form exists 
relating the linear and nonlinear correlation functions. The 
method was extended to power spectra by Peacock & Dodds 
(1994) and modified by Jain et al. (1995) to cope with steep 
power law fluctuation spectra, P(k) oc k n , with n < — 1. 
We have found that the formula of Jain etal. gives more 
self-consistent results for the n ~ — 2 linear power spectra 
discussed here than the revised formula given by Peacock & 
Dodds (1996). 

The linear to nonlinear transformation is given by 

A 2 (k NL )/b(n) = f NL [A 2 L (k L )/b(n)} (1) 
k L = [l + A 2 NL (kN L )]' 1/3 k NL , (2) 

where the subscripts L and NL refer to linear and nonlin- 
ear respectively and A(fc) = 47rfc 3 P(fc)/(27r) 3 is the frac- 
tional variance of the density field in bins of In k. The factor 
b(n) = [(3 + n)/3] l :3 is a function of the effective spectral 
index of the density fluctuations, defined as the local slope 
of the linear power spectrum at the scale on which the vari- 
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^Missing because of space restrictions in astro-ph. Find it in: http://star-www.dur.ac.uk~cmb/project.htm] 



Figure 2. A comparison of the evolved density field in the APM2(a) (top) and CDM2(a) (bottom) simulations, which were started with 
the same random phases. The density is binned on a 256 3 grid and is smoothed with a Gaussian filter to blur the pixels. The greyscale 
shows the logarithm of the density. The slices are 3h~ Mpc thick and 600?i — 1 Mpc square. 
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ance is unity. Using the functional form for the inverse of 
/jvl given by equation 7(b) of Jain el al. (1995), the lin- 
ear power spectrum corresponding to the measured APM 
galaxy power spectrum can be calculated iteratively. The 
linear APM power spectrum is shown by the solid line in 
Figure jj], with the measured APM galaxy power spectrum 
shown by the open circles. The errorbars show the la scat- 
ter in the mean from averaging over the APM Survey spilt 
up into four zones (BE93, BE94). The linear APM power 
spectrum is smoother than the measured spectrum and is 
better fitted by a simple analytic form; for k < 0.6/iMpc - 

PAP M (k) oc — , (3) 

[l+(fc/fc c ) 2 ] 3/2 

with k c « 150 H /c (BG96). 

The linear APM power spectrum is used to generate the 
initial density fluctuations in a N-body simulation. The sim- 
ulation is evolved until the variance measured in spheres of 
radius 30/i _1 Mpc matches that in the APM Survey. Several 
sets of simulations with APM initial conditions are used in 
this paper. APM1 consists of one simulation with 160 3 par- 
ticles in a 440/i _1 Mpc box. APM2 has one realization with 
200 3 particles in a 600fo _1 Mpc box and APM3 is an ensem- 
ble of five simulations (a)-(e), with half as many particles 
as APM1 and a slightly smaller box. The parameters of the 
simulations are listed in Table [j]. The power spectrum of 
the evolved simulation APM2(a) is shown by the dashed 
line in Figure [l]. The evolved power spectra give a very close 
match to the measured APM power spectrum. In all cases, 
we generate the initial conditions using a FFT on a 256 3 po- 
tential grid (N g = 256). The softening length of the APM3 
and CDM3(a) runs was adjusted to be comparable to that 
used in the APM1 run. All simulations were run using the 
P 3 M particle-particle/particle-mesh code of Efstathiou et al. 
(1985). 

We have run two simulations with the same random 
phases as APM2(a) and APM3(a), but with the standard 
CDM power spectrum, Q = 1 and h = 0.5; these are 
called SCDM2(a) and SCDM3(a). We also ran another CDM 
simulation with the same random phases as APM3(a) and 
SCDM3(a), but with a low density parameter and a nonzero 
value of the cosmological constant, f2 = 0.2, A = 0.8 and 
h — 1, which is called LCDM3(a). The initial density field 
in the CDM simulations is set up using the transfer func- 
tion of Bond & Efstathiou (1984) for a universe with baryon 
density Qb = 0.03. This transfer function can be expressed 
in terms of a parameter P = fi/i (Efstathiou, Bond & White 
1992); note that this definition of the shape parameter V is 
relative to a model with f2s = 0.03 and differs slighty from 
that adopted by Peacock & Dodds (1994). In all cases we 
have run the CDM simulations so that the linear variance 
on scales of 8/i _1 Mpc is as — 0.84 (note that this value 
does not take into account any evolution in the clustering, 
and corresponds to clustering at the mean redshift in the 
APM, e.g. G95). The SCDM simulation has more power on 
small scales and less power on large scales than the APM 
run. This can be seen in a comparison of the particle dis- 
tributions from APM2(a) and SCDM2(a) shown in Figure 
|i| The figure shows a slice from the simulation box, after 
the particle density has been tabulated on a 256 3 grid and 
smoothed on small scales with a Gaussian filter. The slice 
shown is ~ 3/i -1 Mpc thick and 600ft _1 Mpc square. 
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Figure 3. Comparison of the theoretical (smooth curve) and 
measured counts (histogram) in radial shells for two mock cata- 
logues made from the simulations SCDM3(a) (top) and APM3(a) 
(bottom). 

3 MOCK APM MAPS 

We transform the N-body simulation into a mock APM cat- 
alogue of angular positions by the following steps: 

(i) Select an arbitrary point in the simulated box to be the 
local 'observer'. 

(ii) Apply the APM Survey angular mask, including plate 
shapes and holes. 

(iii) Include a simulated particle at coordinate distance x 
from the observer with probability given by the selection 
function ip{x). 

The discreteness of the density field in the N-body sim- 
ulations means that the final maps have a slightly lower 
density than the real APM map. The total number of parti- 
cles is about 8 x 10 5 compared with 1.3 x 10 6 galaxies in the 
APM Survey to the same apparent magnitude limit. This 
introduces additional shot-noise in the measurements which 
is corrected in the standard way (e.g. G94). The simulations 
use a periodic box, so we replicate the box to cover the to- 
tal extent of the APM volume (over 1200/i -1 Mpc, beyond 
were the expected number of galaxies is of order unity). By 
comparing the results from different box sizes we have veri- 
fied that this replication of the box does not introduce any 
spurious correlations on large scales. 

3.1 The selection function 

The selection function ip(x) is the normalized probability 
that a galaxy at coordinate distance x is included in the 
catalogue. This probability is proportional to the estimated 
number of galaxies at this coordinate: 

i>{x) = r dq 4>{q) (4) 

Jgi(x) 

where ip* is adjusted so that the probability integrates to 
unity over the sample. (f>(q) is the luminosity function and 
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qi{x) and qz(x) are the scaled luminosities corresponding to 
the lower and upper limits in the range of apparent mag- 
nitudes used to build the galaxy sample or catalog under 
study. In our case these are bj = 17 and bj = 20 respec- 
tively. G95 constructed a x 2 test to find contours of the 
values of the luminosity function parameters that best fit 
observational constraints on the luminosity and redshift dis- 
tribution; the redshift evolution of the luminosity function as 
parameterised as <j>* = 4>o(l + <fi*z); a = ao + a.\z and M = 
Mq + M^z. Here we use the best fit parameters obtained 
by G95: <t>\ ~ 0, ax = -4 and Ml = -2 and the zeroth- 
order values of Loveday etal (1992): 4>* Q = 0.0112/i 3 Mpc~ 3 , 
Mq = —19.73, ao = —1.11. BE93 proposed a functional 
form for the redshift distribution N(z), discussed below in 
Section ^. This N(z) distribution gives very similar results 
for the selection function. 

Figure ^| shows a comparison between the expected 
number of galaxies, n(y)y 2 Ay, at different radial depths (in 
comoving coordinates y) given by the input selection func- 
tion compared to the measured counts for two different mock 
catalogues. 



3.2 Equal area projection maps 

We have made equal area projected maps from the mock cat- 
alogues. To facilitate a comparison between the maps made 
from the different simulations and the map of the APM 
Galaxy Survey, the maps have been turned into greyscale 
plots shown in Figure ^. Grey intensity increases as a low 
power (~ 0.1) of the point density. The mock catalogues are 
from the same realization of the random seeds, and there- 
fore have the same fluctuations in the same places but with 
different amplitudes, given by the difference in the initial 
power spectrum and its subsequent non-linear evolution. In 
the notation of Table |^ these maps are from the simulations 
APM3(a), SCDM3(a) and LCDM3(a). The real APM map 
has been diluted to show the same mean surface density. 
The angular correlations are given in Figures |H] and ^. A 
visual comparison shows that the SCDM model does not 
have as strong large scale fluctuations as the APM map, 
which is confirmed by Figure ^ (as found earlier by Maddox 
et al. 1990c). The SCDM distribution is quite smooth on 
the largest scales. One can also see how both CDM mod- 
els have larger fluctuations on the smallest scales in these 
maps, showing a distinctive granulation in grey scale. The 
mock APM map is the closest of the models to the real cat- 
alogue, as expected from the very good agreement in the 
variance (c/Figure Of course, the locations of individual 
structures in the real and mock APM maps do not coincide. 
Any statistical differences are due to differences in higher 
order correlations. 



3.3 Angular correlations 

Figure ^ compares the variance W2 =< 8 2 > c , of angu- 
lar fluctuations S in cells of radius 9. The angular variance 
in the APM Survey is shown by the points with errorbars 
(G94). The lines show the variance in the mock catalogues 
made from the SCDM3(a) (shot-dashed), LCDM3(a) (long- 
dashed) and APM3(a) (solid) simulations. The CDM mock 
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Figure 5. Comparison of variance of angular counts-in-cells, 
u>2(S), in the simulated mock catalogues with the results in the 
real APM data (symbols with errorbars). A single realization of 
the maps made from the SCDM3(a), LCDM3(a) and APM3(a) 
simulations, normalised to a linear variance of as = 0.84, are 
shown as short-dashed, long-dashed and continuous lines. 
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Figure 6. Comparison of the skewness 1D3 of angular counts- 
in-cells, in the simulated mock catalogues with the results in 
the APM Survey (symbols with errorbars). The SCDM3(a), 
LCDM3(a) and APM3(a) simulations, normalised to linear as = 
0.84, are shown as short-dashed, long-dashed and continuous 
lines. 
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Figure 4. Comparison of equal area projections of maps made from simulated catalogues with the real APM Galaxy Catalogue (Maddox 
et al. 1990a, b,c) (top). The surface density of galaxies is represented by a greyscale, with the densest regions being the brightest. In each 
map, the same total number of galaxies and the same greyscale calibration are used. The maps extend about 120 degrees in RA and 60 
degrees in DEC, covering about 20% of the southern galactic cap, with a mean depth of 400/i _1 Mpc. The 185 overlaping square UK 
Schmidt plates in each map correspond roughly to 5 degrees on a side. All maps have similar amplitudes of fluctuations (1S2) at 1 degree. 
^From top to bottom we show the real APM Survey, the standard CDM map made from SCDM3(a), a lambda-CDM map made from 
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Figure 7. Comparison of variance of counts-in-cells, £ 2 (ii!), i n 
the simulated box (lines) with the results inverted from the cor- 
responding angular mock catalogue (points) . The right panel cor- 
responds to the T = 0.5 CDM model SCDM3(a), while the left 
panel shows the results from the APM3(a) mock catalogue. 



catalogues all have, by construction, the same linear the- 
ory normalisation as ~ 0.84, but the non-linear as is, of 
course, slightly different in each case. The simulated APM 
mock catalogue is slightly off the measured APM values 
around 2 — 4 degrees, this is also true for the mean of differ- 
ent realizations, and is probably due to slight inaccuracies 
in the match of the evolved power spectrum in the simu- 
lation to the measured galaxy power spectrum. Note that 
all catalogues have similar power at around 1 deg. As ex- 
pected the LCDM model matches well the shape of the vari- 
ance at larger scales while the SCDM does not have enough 
large scale power. At smaller scales the CDM models have 
too much power, as noted previously (e.g. Efstathiou et al. 
1990b, G95, BG96). 

Figure |(| shows the skewness or third order reduced 
cumulant H13 =< S 3 > c , for the same single realization of 
each model compared to the APM measurements. The errors 
show that at scales bigger than 1 degree there are very large 
sampling fluctuations in W3. This is more dramatic in single 
realizations of each mock catalogues which cover a smaller 
volume than the real APM. It is therefore dangerous to draw 
any conclusions, from this figure alone, at scales 8 > 1 deg. 
By comparing different realizations, we note that the mean 
at 8 > 1 deg. comes closer to the APM observations. 

In the following two Sections, we first give a brief review 
of the deprojection algorithms, we then estimate the three 
dimensional statistics from the full simulation box i.e. the £ 7 
using the counts in cells method for the whole simulation box 
(as in Baugh, Gaztanaga & Efstathiou 1995) or the power 
spectrum. The two dimensional measurements of clustering, 
the wj or the angular correlation function, are estimated 
from the mock catalogues (as in G94) and the deprojection 
algorithms described in the previous Section are applied. 



4 RECOVERY OF THE MOMENTS OF 
COUNTS IN CELLS. 

Here we use a simple method for recovering the 3-D variance, 
£ 2 (i?), and higher order reduced moments, £j(R), from the 

2- D correlations, wj(8). This method was introduced and 
applied to the APM Galaxy survey in G94, G95, where a 
full description can be found. 

In a scale-invariant model £ 2 ^ R ' with slope 7, we 
can use the expressions in G95 to relate the estimated an- 
gular amplitudes to the underlying three dimensional am- 
plitudes, i.e. al = £ 2 (7? = 8) and Sj = £j/£ 2 Here we 
consider a distribution that is not exactly scale-invariant 
but has a slope 7 which is a slowly varying function of scale. 
We call this a quasi-scale-invariant model (see G95). It is 
then possible to apply a local inversion at each scale. In 
principle the correlations on all scales R contribute to the 
correlations on angular scale 8, but because the sample has 
a finite depth, T>, there is a characteristic scale R ~ T>8. In 
our analysis we relate angular scales 8 to 3-D scales using 
R — T>8, where T> is the estimated distance which corre- 
sponds to the mean redshift of the sample (see also Peebles 
1980). Although there is some ambiguity as to what the best 
definition of 2? should be, in the scale-invariant regime, we 
find that the estimated amplitudes of £j are insensitive to 
changes in our chosen value of D . 

Thus at each given given scale 8 with local slope 7 = 
7(6), we use the scale invariant expressions to relate the 
estimated local angular amplitudes to the underlying three 
dimensional values. This results in an estimation for ^ as a 
function of the scale R — T>8. This model was used in G94 
and G95 to recover the 3D correlations in the APM Survey. 

4.1 Test of the variance 

Figure [?] shows the inversion of £, 2 (R) from a standard 
r = Qh = 0.5 CDM mock angular catalogue (right panel), 
and for a mock APM catalogue (left panel) compared to 
the corresponding variance £ 2 (-R) estimated directly in the 

3- dimensional simulation [e.g. SCDM3(a) and APM3(a)]. 
The variance recovered from the angular distribution is a 
very good match to the variance measured from the full 
simulation. There is a slight disagreement at scales around 
R ~ 20/i~ 1 Mpc, where there is a rapid change in the slope, 
as expected, but the discrepancies are within la. 

4.2 Test of higher order moments 

The simulations we use have values of Sj which show a small 
variation with scale, e.g. 5*3 oc R~ a , with a ~ 0.1. This 
indicates that strictly speaking neither the scale-invariant 
nor the quasi-scale-invariant models should be used, as Sj 
should be constants in the hierarchical model. Nevertheless, 
we still find reasonable agreement from the inversion when 
we compare local values of Sj. 

Figures ^ and ^ show the inversion of Ss(R) and Si(R) 
from the standard T = 0.5 CDM mock angular catalog and 
for the APM-like mock catalogue compared to the corre- 
sponding amplitudes estimated directly in the 3-dimensional 
simulated box. At scales 20/i _1 Mpc > R > 6/i _1 Mpc, the 
amplitudes recovered from the angular distribution are in 
good agreement with the original amplitudes. At larger 
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Figure 8. Comparison of the skewness, S%(R), in the simulated 
box (lines) with the results inverted from the corresponding an- 
gular mock catalogue (points). The top panel shows the T = 0.5 
SCDM3(a) model, while the bottom panel corresponds to the 
APM3(a) simulation. 
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Figure 9. Comparison of the kurtosis, S4,{R), in the simulated 
boxes (lines) with the results obtained from the angular cata- 
logues (as in Figure ph. 



scales, sampling fluctuations are very large, whereas at 
smaller scales, there are some systematic differences which 
seem more important for the APM model, which has a 
steeper power spectrum. As expected, the inversion method 
seems to work better for distributions where the Sj are 
closer to being constants, e.g. SCDM. Note that the mea- 
sured amplitudes Sj in the APM are closer to a constant 
than either of the models we study here (G94,G95) and one 
would then expect an even better agreement in this case. 
The discrepancies at small scales could also be due in part 
to shot-noise in either the angular or the 3-dimensional dis- 
tribution. 

As pointed out in Gaztanaga & Bernardeau (1997), 
there are several effects that make this type of comparison 
difficult. First, volume and boundary effects are important 
on scales > 2 deg and tend to produce smaller values of the 
projected amplitudes S3 and S4. Second, the simple hierar- 
chical model for projections commonly used in the literature 
(e.g. by Groth & Peebles 1977, Fry & Peebles 1978) is not 
accurate on quasi-linear scales, as indicated in Bernardeau 
(1995). These two effects compete with each other and it is 
not clear how the projection model should be improved to 
allow a better reconstruction. 



5 RECOVERY OF THE POWER SPECTRUM 

BE93, BE94 developed an iterative technique to numerically 
invert Limber's (1954) equation which relates a measure of 
clustering in 2D to an integral of the 3D power spectrum 
multiplied by the survey selection function. BE93 used the 
measured angular correlation function of the APM Survey, 
w{0) to obtain an estimate of the 3D power spectrum, whilst 
the 2D power spectrum, P2{k) was used in BE94. An esti- 
mate of the real space correlation function has also been 
made in the same way (Baugh 1996). 

This algorithm for the numerical inversion of Limber's 
equation does not rely upon the initial form chosen for the 
power spectrum and can reveal features that would be dif- 
ficult to parameterize in a simple way. The technique is nu- 
merically stable, unlike the use of Mellin transforms which 
involve differentiation of noisy quantities (Fall & Tremaine 
1977), and it has been shown to rapidly converge to stable 
solutions (BE93). 

The integral equation relating w(9) to the 3D power 
spectrum, P(k), is given by (BE93, see Peacock 1991 for the 
non-relativistic form) 



w(u>) 



P(k)kg(kuj)Ak, 



(5) 



where the angular variable is u> = 2sin(#/2) and the kernel 
function is an integral over the survey selection function 



%£J (kux)dz, 

where F(x) depends upon the cosmological model (see Pee- 
bles 1980, §56) and Qs is the solid angle of the survey. The 
time evolution of the power spectrum is parameterised as 
P(k,z) = P(k)/(1 + z) a , where a = corresponds to the 
pattern of clustering being fixed in comoving coordinates, 
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which is the case we use in this paper. This is a necessary 
oversimplification as we have an observed quantity that is 
a function of only one variable. Furthermore, the median 
redshift of the APM Galaxy survey is z m ~ 0.12 and the 
corrections for redshift evolution are small. 

The redshift distribution of survey galaxies is parame- 
terised as (BE93) : 



3JV(m)n s 2 
2zl 2 



/ 2 \ 3/2 

exp(-(-) )d 



with the median redshift given by 
z m = 1.412.!s c = 0.016(foj 



17) 1 ' 5 +0.046 



(7)t 



(8); 



for apparent magnitudes b.j > 17. This form was chosen to 
provide a fit to the redshift distribution in the Stromlo/APM 
survey (Loveday et al. 1992) and to the fainter surveys of 
Broadhurst et al. (1988) and of Colless et al. (1990, 1993). 
Redshifts have now been measured for galaxies in the mag- 
nitude range covered by the APM Survey, 17 < bj < 20, 
(Ellis et al 1996) and the redshift distribution is in good 
agreement with the form that we have adopted (Efstathiou 
private communication) . 

The r th iteration of the Lucy algorithm gives an esti- 
mate of the data, of 



u r (uji) = S P r (kj)g{kjU>j)kj Alnfc 



(9) 



which is compared with the 'true' data, w°(6) in order to 
generate a new estimate of the power spectrum: 



^gikjUi) Ainu 



(10) 



The summations have typically 60 logarithmic bins for the 
data and 30 logarithmic bins for P(k) in the range 3xl0~ 3 < 
k < 30/iMpc" 1 . 



5.1 Test of the Recovery of P(k). 

In all cases, unless otherwise stated, we use the mean an- 
gular 2-point correlation function and its variance in 4 in- 
dividual disjoint zones (shown in Figure 2 of BE94) to re- 
cover the power spectrum. The results of the inversion of 
equation ^ are illustrated in Figure [H^, for the two CDM 
models SCDM3(a) and LCDM3(a), which have power spec- 
tra with very different amplitudes and curvatures at a given 
wavenumber. There is very good agreement in each case, up 
to the largest wavenumbers sampled in the simulation box in 
these runs, k ~ 2tt/L ~ 0.015ftMpc _1 . In section §[u| below 
we present a more detailed comparison for larger scales. 
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Figure 10. Symbols with errorbars show the recovered 3D P(k) 
from the inversion of equation ^ for maps made from the CDM 
simulations: a) LCDM3(a) b) SCDM3(a). The dotted and dashed 
lines show for reference the linear power spectra of the LCDM and 
SCDM models respectively. The solid line in each case shows the 
nonlinear P(k) measured from the simulation. 



6 MEASURING THE BREAK IN P(K). 
6.1 Accuracy on large scales 

In order to show that the inversion method can accurately 
recover features at small wavenumbers (large scales), such as 
the break in P(k), we now concentrate on the largest volume 
simulations, with a box size of L = 600ft _1 Mpc (see Table 
^|). We study a single mock angular map from SCDM2(a) 
and APM2(a), with the same phase correlations and position 
for the observer. Figure [ll] shows a comparison of the initial 
linear P(k) (dashed lines) with the non-linear P(k) in 3D 
from the full box (continuous line) for both: (a) SCDM (right 
panel) and (b) the APM model (left panel). Note how even 
the P(k) measured in the 3D box has large fluctuations at 
small k, and in particular a large spike at k ~ 0.03ftMpc _1 . 
This is due to the small number of modes available to esti- 
mate P(k) on these scales with the Fast Fourier Transform 
(FFT) technique. These estimates have not been averaged in 
bins and the initial spectrum amplitudes are drawn from a 
Gaussian distribution (and are not set equal to the mean). 
The mode where this spike is located only corresponds to 
nx=2 ny=2 nz=0, so that there are few modes to average 
over. This is seen in both the APM and CDM P(k) in this 
plot, due to these simulations being set up with the same 
phase distributions. 

On large scales in Figure |ll| we plot the power spectrum 
at the individual Fourier modes. At large wavenumbers we 
have binned the 3D FFT estimation for clarity. The recov- 
ered P(k) from the angular two-point function (points with 
errorbars) shows excellent agreement with the original P(k). 
Hence the volume of a sing le N-body box (L = 600ft -1 Mpc) 
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Figure 11. Comparison of the recovered 3-D P(k) (points with 
errorbars), from the angular catalogues, with the P(k) from the 
corresponding 3D simulation (continuous line). The left hand 
panel shows the results for APM2(a) and the right hand panel 
shows SCDM2(a). The long and short dashed lines correspond to 
the r = 0.5 linear CDM P(k) and the APM linear power spec- 
trum respectively. 




0.01 0.1 1 0.01 0.1 1 
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Figure 12. Comparison of the recovered 3-D local slope n(k) 
(points with errorbars), from the angular catalogues, with the 
slope in the corresponding 3D power spectrum (continuous line). 
The long and short dashed lines correspond to the T = 0.5 linear 
CDM model and the APM linear model. The left and right panels 
show the results in the small (400?t — 1 Mpc) and large (600/i _1 Mpc 
boxes). The top and bottom panels correspond to SCDM and 
APM simulations respectively. 



is large enough to simulate and recover large scales features 
in P(k), even at k ~ 0.01/iMpc" 1 . 

It is clear from this figure alone that there is a significant 
measurement of the break of the power spectrum. To make 
this more qualitative we now turn to the local slope of P(k). 

We want to focus in more detail on the shape of the 
power spectrum by estimating the local logarithmic slope: 



n(k) 



dlogP(fc) 
d log k 



(11) 



To do a numerical estimation we first bin the P(k) data 
and use standard polynomial interpolation and numerical 
differentiation (e.g. Press et al. 1992) in logarithmic space. 
The error in the slope is obtained assuming no spread in k: 



An(fc) 



dAP(k)/P(k) 
dlog k 



(12) 



This approach seems to work well in the mock maps and 
avoids spreading the systematic errors coming from the sam- 
pling variance, which typically introduces a larger uncer- 
tainty in the amplitude of the correlations than in their 
shape (see Figure 4 in Baugh etal. 1995). 

Figure jl2] shows the results for single realizations of 
the SCDM and APM models for two different box sizes; 
SCDM2(a) and AMP2(a) with box size L = 600/i _1 Mpc 
and SCDM3(a) and APM3(a) which have a box size of 
L — 400/i~ 1 Mpc. The largest scales sampled in each pair 
of simulations correspond to wavenumbers of k ~ 2tt/L ~ 
0.01/iMpc" 1 and k ~ 0.015/iMpc" 1 respectively for the 
L = 600/i _1 Mpc and L = 400/i _1 Mpc boxes. The small- 
est scales sampled are limited by the Nyquist frequency of 



the FFT grid (of size N g ): k ~ N g n/L, or for large enough 
N g by the numerical resolution (e in Table |l]). 

Figure ^ shows that the recovered slope matches 
closely that obtained directly in three dimensions, both 
in the non-linear (k > 0.2foMpc -1 ) and linear (k < 
0.1/iMpc -1 ) regimes. The particular realisations of the 
smaller boxes shown in the Figure [l2] have a flatter slope 
on large scales in three dimensions than the corresponding 
linear spectrum due to finite volume effects. This effect is 
also reproduced in the recovered slopes. 

The break in P(k) corresponds to n = 0, and is well 
traced within the errors in the larger boxes. 



6.2 Implications for the APM Power 
Spectrum. 

Table |2| shows the values of P(k) recovered from the two- 
point correlations in the APM angular Galaxy Catalogue. 
These are essentially the same as in Figure 7 of BE93, al- 
though there are small differences corresponding to a dif- 
ferent number of iterations in the Lucy algorithm (chosen 
here to provide the minimum \ 2 match to the angular cor- 
relation function) . Figures [l^ and jl4] illustrate the implica- 
tions of our findings for the power spectrum recovered from 
the APM. Figure [13 shows the reconstructed slope in the 
APM Galaxy power spectrum, while Figure jl4] shows the 
corresponding P(k). Symbols with errorbars correspond to 
the mean and variance in 4 individual disjoint zones (shown 
in Table ^). The break at n — is found to lie between 
k = 0.02 — 0.06/iMpc -1 (between the vertical dotted lines 
in the Figure). This can also be shown directly in Figure 
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Figure 13. The local slope of the power spectrum estimated from 
the APM catalogue. Symbols with errorbars correspond to the 
slope estimated in 4 individual disjoint zones (subsamples). The 
continuous line is from an inversion of the angular correlation 
function measured from the full APM map. The short dashed 
line corresponds to the inversion after subtracting 10~ 3 from the 
angular correlation function in the full APM map. The two long- 
dashed lines correspond to linear CDM models with Q.h = 0.5 
(top) and Qh = 0.2 (botom). 



where P(k) shows a significant break on similar scales 
(also bounded by dotted lines). Note that the errorbars are 
comparable in the mock and the real catalogues. 

The power spectrum recovered from the mock cata- 
logues agrees well with the P(k) measured from the unpro- 
tected simulation box, indicating that the volume of a single 
box (L = 600ft" 1 Mpc) is large enough to realize and recover 
a break on scales around k ~ 0.05/tMpc -1 , without any fi- 
nite volume effects. Thus the volume traced by the APM 
Survey (which extends radially well beyond 600ft _1 Mpc) is 
large enough to allow a measurement of the break in the 
power spectrum, n — 0. 

In an extensive analysis of the systematic errors in- 
volved in plate matching, Maddox et al (1996) have placed 
an upper limit of 5w{6) ~ 1 x 10~ 3 on the likely contri- 
bution of the systematic errors to the angular correlations. 
In Figures and[y| the inversion result using the angular 
correlation function measured from the full survey is shown 
as a continuous line. The short dashed lines in these fig- 
ures show how this result for the power spectrum changes 
when an offset of 10 -3 is subtracted from the angular cor- 
relation function in the full APM map. In principle, results 
from individual zones (symbols with errors in the Figures) 
could be affected more by the zone boundary than results 
from the full survey, though on the other hand large scale 
noise from plate matching could be more important for the 
whole survey than for individual zones. For the mock cat- 
alogues, the smaller size of individual zones does not seem 
to introduce important errors at the scales under consider- 
ation (e.g. Figure [li]). Thus, while there is no clear reason 
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Figure 14. Comparison of the recovered P(k) (points with er- 
rorbars), from the angular catalogues with the linear APM model 
(long-dashed line). Panel (a) corresponds to a mock APM cat- 
alogue [APM2(a)]. Panel (b) shows the estimated APM P(k) 
from measurements in the real galaxy catalogue. In both cases 
the points and errorbars correspond to the mean and variance in 
4 individual disjoint zones (subsamples). The continuous line in 
panel (b) corresponds to the inversion result obtained using the 
angular correlation function measured from the full APM map. 
The short-dashed line corresponds to the inversion result after 
subtracting an offset of 10 — 3 from w(9). 



to prefer the estimate of the power spectrum made from the 
full survey to that made from the zones, the later is less 
likely to be affected by any large scale plate matching er- 
rors. As the variance from the different zones includes all 
the above sources of potential error, we take this estimation 
and variance as our best mean and errors. 

Figure [iB] shows the effects of nonlinear evolution in 
the mass power spectrum for the F = 0.5 standard CDM 
model and for two variants of SCDM with T = 0.2. Again, 
we use the form of the CDM power spectrum given by Bond 
& Efstathiou (1984), which is valid for a universe with a 
small baryon density, Qb = 0.03, and we follow the defini- 
tion of T = tth = 0.5 for SCDM adopted by Efstathiou etal 
(1992). For I = 0.5 we show linear theory power spectra 
(solid lines) for two different normalisations to the variance 
in spheres of radius 8/i _1 Mpc; the amplitude of tempera- 
ture fluctuations in the microwave background gives a value 
as ~ 1.2 (e.g. Stompor etal 1995, Bunn, Liddle & White 
1996), whilst normalisation to reproduce the abundance of 
rich clusters requires as ~ 0.50, virtually independent of the 
shape of the power spectrum for f2 = 1 (Eke, Cole & Frenk 
1996; White, Efstathiou, & Frenk 1993). The dashed lines 
give the corresponding predictions for the nonlinear spec- 
tra, using the transformation of Peacock & Dodds (1996) 
rather than Jain et al. (1995), which is not so accurate for 
CDM models (see BG96). The lower set of curves in Fig 
|lf](c) show a critical density model with a Hubble constant 
H = SOkms^Mpc" 1 , but with the SCDM transfer func- 
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Figure 15. The effects of nonlinear evolution on the shape of CDM power spectrum. The solid lines show the linear theory power 
spectrum: the lower and upper curves in (a) and (b) arc for normalisations of as = 0.50 and as = 1.21 respectively in the T = 0.5 
CDM model. The lower solid curve in panel (c) shows a T = 0.2 CDM model for a universe with f2 = 1 and a Hubble constant of 
Ho = 50kms~ 1 Mpc~ 1 normalized to match the COBE result with ag = 0.42. The upper solid curve in panel (c) corresponds to a open 
universe with Q = 0.2 and Ho = 100kms~ 1 Mpc~ 1 , with a normalisation that reproduces the abundance of rich clusters, as = 1.07. 
The dashed curves show the corresponding nonlinear power spectra; in this case we have used the transformation of Peacock and Dodds 
(1996), which is expressed in a form that makes it readily appliciablc to different cosmological models. 

The APM galaxy power spectrum is shown by points with error bars, where we have divided by a bias parameter, b s , squared (values 
indicated in the figure) to match the amplitude of the mass spectrum at different scales: (a) the bias parameters have been chosen to 
match the amplitude of the mass T = 0.5 CDM P(k) on small scales k > 0.5; (b) the match is made to the amplitude of the mass T = 0.5 
CDM spectrum at the scale of the break in the APM power spectrum, k ~ 0.05; (c) match to the amplitude of the mass T = 0.2 CDM 
spectrum on small scales, k > 0.5.. At small wavenumbers, k < O.Ol/iMpc -1 , the estimate of the APM power spectrum is dominated by 
systematic and random errors in the catalogue 



tion altered by using T — 0.2. The normalisation of this 
curve matches the COBE detection, with as = 0.42. The 
upper curves in (c) are for an open model with density pa- 
rameter O = 0.2 and Ho = 100kms _1 Mpc _1 . In this case, 
the model is normalised to reproduce the abundance of rich 
clusters with cr 8 = 1.07 (Eke etal 1996). 

The APM galaxy power spectrum has been plotted to 
match the amplitude of the mass power spectrum, divid- 
ing the observed amplitude by a bias parameter squared. 
In the (a) panel, the bias parameters have been chosen to 
match the galaxy power spectrum to the amplitude of the 
mass power spectrum at wavenumbers 0.3 < k < lOftMpc -1 , 
whilst in the (b) panel, the match is made at the scale of the 
break in the APM power spectrum. Panel (c) shows the case 



for r = 0.2 normalized to small scales. This Figure demon- 
strates the basic problem of the CDM model; the shape of 
the power spectrum cannot be made to match to observed 
galaxy power spectrum at both large and small scales, unless 
some complicated biasing prescription is invoked, in which 
the bias would need to vary significantly with scale. 

6.2.1 The effect of biasing 

The fluctuations traced by the galaxy distribution might be 
different, or biased, from the underlying mass fluctuations 
(e.g. Bardeen et al 1986). We will argue here that the effect 
of this biasing is not important for the shape of the APM 
power spectrum at large scales. 
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Assume that the (smoothed) galaxy fluctuations S g are 
related to the mass S m fluctuations by a local transforma- 
tion: 5 g (x) — F[S m (x)], and that this relation can be given 
as a Taylor series: F = bi5 m -\-b28m + ... . Then the two-point 
function £f (r) =< 5 g (x)S g (x + r) > on a scale r will be just 
given by: 

&(r) = b{&(r) + &1&2 < * m ("0M* + rf > + 

+6162 < 5 m (x)S m (x + r) > +■■■ (13) 

were all further terms are of order 4 or greater in 5 m , 
and therefore correspond to either higher order correlations, 
£j with J > 2, or higher powers in £2. If S m is Gaussian or 
hierarchical (as in the case for gravitational evolution) the 
higher order correlations £j are at most of order ■ This 
means that at large scales, where £2 < 1, the first term is 
the dominant one, so that only the amplitude but not the 
shape of the two-point statistics is changed by biasing. This 
effect have been found in N-body simulations and toy biasing 
models (Weinberg 1994, Mo, Jing & White 1997, Gaztanaga 
& Lacey 1997). 

For the APM power spectrum a wavenumber around 
k ~ O.lftMpc" 1 , corresponds to a top-hat radius of 7? ~ 
%/k ~ 30/i _1 Mpc. For any reasonable biasing model relat- 
ing galaxy fluctuations, S g , to the underlying matter fluc- 
tuations, 5 m , the matter density fluctuations are very small 
around R ~ 30/i _1 Mpc. The independent constraints on 
the normalisation of mass fluctuations discussed above give 
values of around unity for the variance in spheres of ra- 
dius 8/i _1 Mpc. To have rms fluctuations of order unity at 
R ~ 30/i _1 Mpc would imply as > 3. 

Thus from the above arguments, the small variance on 
large scales, R > 8/i _1 Mpc, means that it is reasonable to 
assume that the galaxy shape of P(k) for k < 0.1/i.Mpc -1 
corresponds to the shape of the underlying linear matter 
power spectrum. This argument, just based on the small- 
ness of the variance and the hierarchical structure, can also 
be applied to gravity, as the leading contribution to the cor- 
relation functions in perturbation theory is indeed exactly 
given by a local transformation (see Fosalba & Gaztanaga 
1997). This is clearly illustrated in Figure [l2| By comparing 
the linear and non-linear shape of P(k), one case see that 
it has not been changed significantly by gravitational evo- 
lution on scales where the rms fluctuations are small, i.e. 
k < 0.1/iMpc -1 . 

6.2.2 Variations of CDM models 

A simple variation of CDM models is to introduce a tilt in 
the initial power spectrum so that: P(k) — k n °T(k), where 
T(k) is the transfer function (e.g. Bond & Efstathiou 1984, 
Bardeen et al. 1986) and no, is the primordial spectral index, 
no 7^ 1. Unless the tranfer function T(k) somehow depends 
strongly on no, the local slope of a given tilted CDM model is 
similar to that of the corresponding standard scale invariant 
model (where no = 1), given by n = no + dlog(T) / dlog(k) , 
with the shift due to the tilted value of no. Thus, tilted mod- 
els can only scale up or down the CDM predictions in Figure 
|l3| , and therefore can not account for the APM observations. 

The measurement of the abundance of deuterium in 
high redshift hydrogen clouds is provoking much debate in 



Table 2. Values of the estimated power spectrum P(k) recovered 
from measurements in the (real) APM angular galaxy catalogue, 
corresponding to the mean and error from the variance in 4 indi- 
vidual disjoint regions in the catalogue. 



k 


P(k) 


AP(fc) 
(h- 3 Mpc 3 ) 


/iMpc" 1 


(h- 3 Mpc 3 ) 


0.0032 


7198 


4345 


0.0043 


6891 


3278 


0.0060 


5805 


2126 


0.0082 


5386 


1543 


0.0113 


6158 


1620 


0.0155 


8134 


2026 


0.0213 


10174 


2803 


0.0292 


10251 


3682 


0.0401 


9821 


3232 


0.0551 


10776 


523 


0.0757 


9440 


1770 


0.104 


6299 


1383 


0.143 


3358 


590 


0.196 


1754 


173 


0.270 


1048 


58 


0.371 


675 


51 


0.509 


451 


41 


0.700 


309 


33 


0.961 


214 


24 


1.32 


146 


16 


1.81 


96.7 


8.8 


2.49 


61.4 


4.6 


3.42 


38.1 


2.4 


4.70 


23.8 


1.3 


6.46 


15.1 


0.7 


8.88 


9.63 


0.43 


12.2 


6.19 


0.25 


16.8 


4.15 


0.15 


23.0 


3.05 


0.10 


31.6 


2.42 


0.07 



the literature (e.g. Rugers & Hogan 1996, Tytler, Fan & 
Buries 1996). Consequently the baryon density of the uni- 
verse is uncertain and possible values fall in a wider range 
than was previously accepted. In the limit of a high baryon 
density (i.e. Qb ~ 0.1), the power spectrum of the mass 
is modified. A full calculation of the transfer function (e.g. 
Slejak & Zaldarriaga 1996) indicates that the high baryon 
density introduces features or 'wiggles' into the shape of the 
power spectrum on large scales (see also Goldberg & Hamil- 
ton 1997 for a discussion of how these peaks could be used 
to constrain the value of SIb). We have used the CMBfast 
code of Slejak & Zaldarriaga to compute the shape of the 
power spectrum in a CDM universe with Qb =0.1 and 
fl = 1. The resulting modification of the power spectrum 
compared with the Bond & Efstathiou (1984) transfer func- 
tion for ils = 0.03 is insufficient to improve the agreement 
with the APM power spectrum. 
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7 CONCLUSIONS 

The algorithms tested here successfully recover the power 
spectrum and higher order cumulants in three dimensions. 
There are no systematic shifts or biases in the inferred cor- 
relations resulting either from the deprojection techniques 
or from the process of projecting the original particle distri- 
bution. The 3D variance recovered from angular catalogues 
is in good agreement with the input model, confirming the 
results in G95. For higher order correlations the deprojec- 
tion method studied here, and also used in G94, Gaztanaga 
& Frieman (1994), G95 and BG96, seems to be adequate, at 
least for intermediate scales, 20/i~ 1 Mpc > R > 6ft~ 1 Mpc, 
although one would in principle expect deviations from 
the smiple hierarchical, according to perturbation theory 
(Bernardeau 1995). A more detailed analysis of this point 
is presented elsewhere (see Gaztanaga & Bernardeau 1997). 

It is possible to recover the detailed shape of the power 
spectrum with errorbars similar to those quoted by BE93. 
As pointed out there (and also in G95), the uncertain- 
ties in the selection function do not have much effect on 
the recovered shape. The deprojection algorithm is able to 
distinguish sharp features, such as the one between k ~ 
0.07 — 0.2/iMpc -1 shown in Figure [l^, first remarked upon 
by BE93. For this range of wavenumbers, the best fitting 
CDM model has flh ~ 0.2, as pointed out by Efstathiou 
etal (1990b) and Peacock & Dodds (1994). However, the 
break in the power spectrum in this particular CDM model 
is broader and at a larger scale than the break in the APM 
power spectrum. 

We have shown that the volume traced by the APM 
Survey is large enough to allow a significant measurement 
of the break in the power spectrum, n = 0, as found on 
scales around k ~ 0.05/iMpc - . We have also shown (See 
Figure 14) that possible systematic errors involved in the 
APM plate matching lie within our estimated errors. 

Peacock & Dodds (1994) report a break in the power 
spectrum at a wavenumber of k ~ 0.03/iMpc -1 using spec- 
tra measured from a range of different surveys. The volumes 
mapped out by these surveys span a considerable range. We 
have found that only our largest simulation boxes allow the 
break to be measured accurately, both in the direct estima- 
tion of the power spectrum in three dimensions and in the 
recovered spectrum obtained from the projected catalogue. 
The size of the largest box we use, L = 600h~ Mpc, is much 
greater than the median depth of any of the redshift surveys 
available to Peacock & Dodds, indicating that finite volume 
effects could have altered the shape of the power spectra 
estimated from individual surveys on large scales (as found 
in Figure |l| for the 400/i _1 Mpc boxes). Other sources of 
uncertainty in this type of compilation include the differ- 
ent selection biases applied, the differences in the intrinsic 
luminosities of the objects selected in the catalogues and 
the large sampling variance from the smaller surveys. Fur- 
thermore, the linearisation process applied to the measured 
power spectra involves a correction for the distortion of the 
pattern of clustering by galaxy peculiar velocities (Kaiser 
1987), which is both model and catalogue dependent (e.g. 
Smith etal 1997). 

The location of the break that we find in the galaxy 
power spectrum matches that found in power spectrum of 
galaxies clusters, both from a compilation based on the Abell 



catalogue (Einasto etal 1997) and from a cafefully selected 
redshift sample drawn from the APM Cluster catalgoue 
(Tadros 1996, Tadros & Dalton 1997, Dalton etal 1992). 
The physical interpretation of the break at 



fc fl ~ 0.05 /iMpc" 1 



150 



Ho 



(14) 



f ound in the APM is unclear. We have argued in section 
§ |6.2.1 that the galaxy shape of P(k) for k < 0.1/iMpc -1 cor- 



responds to the shape in the underlying linear matter power 
spectrum. For inflationary models with Cold Dark Matter 
(CDM) the break in the power spectrum at wavenumber 
ks corresponds to the Hubble radius when the universe be- 
comes matter dominated. This is because the amplitude of 
fluctuations is frozen as they enter the Hubble radius dur- 
ing the radiation dominated era (see Bond & Efstathiou 
1984, Bardeen et al. 1986). The wavelength of the Hubble 
radius at this epoch is As ~ 10(f2/i) -1 /i -1 Mpc (e.g. Kolb & 
Turner 1990), where Q is the total matter density in units 
of the critical density, which corresponds to a wavenum- 
ber of k B ~ 0.1(n/i)/iMpc _1 . Thus, for CDM-like models 
the range of the scales we find for the break in the APM, 
k = 0.02-0.06/iMpc" 1 , implies 0.2 <flh< 0.6. To be more 
precise we perform a x 2 fit to the CDM models in Bond & Ef- 
stathiou (1984) using the four APM P(k) points in the range 
k = 0.02 - 0.06/iMpc" 1 to find Q,h ~ 0.45 ± 0.10 (tth = 0.2 
produces a \ 2 — 9, while fi/i = 0.4 gives \ 2 — 1-3). Thus the 
case fl = 1 requires h ~ 0.45 ± 0.10 while an open universe 
or one with a non-zero cosmological constant, A, can accom- 
modate other values of the Hubble constant h. For purely 
relativistic dark matter, like neutrinos, the scale at which 
the amplitude of fluctuations are damped is typically larger 
than for CDM, corresponding to the Hubble radius when 
the universe becomes non-relativistic. For these models the 
measured break yields correspondingly larger values for Qh. 

As shown in Figure the sharp change in the local 
slope of the APM between k ~ 0.05 - O.lhMpc" 1 is not 
compatible with any CDM model, which have a broader 
peak. Note that in Figure [l^, the results are independent 
of uncertainties in the overall normalization or in any linear 
bias that may be applied, unlike Figure [l^. We have also 
shown that non-linear evolution, is not sufficient to modify 
the shape of the linear CDM power spectrum to provide a 
good match to the shap e of th e observed APM spectrum. 

We have argued in §6.2.2 that simple variation of CDM 
models, such us tilted or higher £Ib models can not account 
for the APM observations. Models in which a large fraction 
of the matter is relativistic (such as Mixed Dark Matter) are 
more likely to match this type of sharp feature. The scale 
found here for the break, around k ~ 0.05/iMpc -1 , could 
give interesting constraints for these models. 
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